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1. Introduction 

Understanding the properties of the vacuum of gauge theories is still an interesting field, po- 
tentially rich of surprises and unexpected applications. 

It is well recognised today that relativistic field models can be used as effective theories to 
describe the low-energy excitations in condensed matter systems; in particular, they can be applied 
to a wide class of planar systems (see for example Ref. part VI). 

Recently much attention has been devoted to the application of the lattice approach to describe 
some properties of a single sheet of carbon atoms arranged in the well known honeycomb structure 
called "graphene" (among the first papers that appeared we want to remind Ref. [||] and Ref. [||]). 

The vacuum structure of lattice gauge theories can be understood probing it by an external 
background field A ext . In this paper, we show that the presence of an external magnetic field 
induces a dynamical symmetry breaking in the particular case of QED in 2 + 1 dimensions with 
dynamical fermions. 

We will also discuss how our results can be applied to describe some striking properties, ob- 
served in experimental studies, of graphene at large external magnetic fields. 

2. The magnetic field background on the lattice 

The introduction of a magnetic background can be done defining on the lattice a gauge- 
invariant effective action r(A ext ) by using the Schrodinger Functional (SF) [H]. The Euclidean 
SF in Yang-Mills theories without matter is defined by 

Z[A^A''] = (A f \e- HT >|A f ) , (2.1) 

i.e. it is the propagation kernel for going from some field configuration A 1 , at time X4 = 0, to some 
other configuration A?, at X4 = T; the operator & projects onto the physical states. 

The lattice SF is given by Z[U^,U l ] = J DUe~ s , where S is the Wilson action modified to 
take in account the boundaries: U (x) X4 =o = U' and U (x) m =t = U*. We define the lattice effective 
action for a background field A ext as 

T(A ext ) = -^ln(^p) , whereZ[[/ ext ]=Z[[/ ext ,[/ ext ]. (2.2) 

Remarkably, it turns out that r(A ext ) is invariant under lattice gauge transformations of the external 
link U ext . Since in this definition = U l , we have periodic conditions in the time direction and, 
due to the lack of free boundaries, the lattice action is now the familiar Wilson action. 
Moreover, it is possible to show that 

r(A ext ) -> £ (A ext ) - £ (0) , (when T -> 00) , (2.3) 

where £"o(A ext ) is the vacuum energy in presence of the external background field. Therefore 
r(A ext ) is the lattice gauge-invariant effective action for the background field A ext . In other words 
to study a theory with an external background field we have to simulate on the lattice the original 
action (i.e. the one without any external field), but introducing proper constraints. 
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In our case we are interested in £/(l) in a uniform external magnetic field H, therefore, after 
imposing spatial and temporal boundary conditions, we have to constrain the spatial lattice links 
belonging to a fixed time slice (for example X4 = 0) to 

Uf xt (x ) = 1 , C/ 2 ext (x) = cos (gHxi ) + i sin (gHxi ) , (2.4) 

where g is the coupling constant. The same constraints are imposed at the spatial boundaries of the 
other time slices, i.e. we require that fluctuations over the background field vanish at infinity. The 
temporal links are not constrained, which is in coherence with the correct definition of the thermal 
partition functional [|p. 

We can see that, since the lattice has the topology of a torus, the magnetic field turns out to be 
quantized: 

2 2k 

agH = —n ext , (next = 0,1,...) , (2.5) 
Lt 

where L, is the temporal lattice extension. A different approach to the problem we have tackled, 
where also a different way to introduce the external magnetic field is implemented, can be found in 
Ref. [§. 

This paper is based on the study of QED with Nf = 2 flavours of 4-component fermions using 
the staggered fermion approach; this means that we need to simulate N = I staggered fermions 
fields x,Z w i tn the Euclidean action [0]: 

N 

S = S g + L YjCi{n)M n ,kXi{k) , (2.6) 

i=l n,k 

and fermion matrix given by 

M n M[U}= £ ^^{[^vWlW-^vWI^J+m^, ( 2 - 7 ) 
v=l,2,3 l 

where tj v (?i) = (— l)"i+ - +"v-i Moreover, we choose the compact formulation of QED, 



S G [U]=P £ l-^(^v(n) + ^ v (n)) 



(2.8) 



where U^ v (n) is the "plaquette variable" and j8 = 1/ (g 2 a), a being the lattice spacing. 

Note that the introduction of the fermions in the theory does not change anything about the 
way we introduce the external field A ext [||]. 

3. Dynamical symmetry breaking 

It is now a general result that a constant magnetic field leads to the generation of a fermion 
dynamical mass []|] in a wide class of 2 + 1 and 3 + 1 dimensional theories; this phenomenon is 
known as magnetic catalysis [JlO|]. 

In the QED in three dimensions it is possible to determine, at first order in perturbation theory, 
the value of the chiral condensate in presence of the magnetic field H (in cgs units) [|ll|] : 

(W = -2Nf\m\c — — V . — . (3.1) 
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After regularisation of the integral and in the limit of small gap Ao = mc 2 , we get: 

ZK ■ J% HceH 

4. Monte Carlo Simulations 

Our simulations have been performed in the weak coupling regime with j8 = 2.0, on volumes 
L 3 , with L = 12, 16,24, for masses in the range 0.005 < rao < 0.05 and for several values of the 
strengths of the magnetic field n ext , namely, n ext = 1,2,3. The Monte Carlo simulation code is 
based on the Hybrid Monte Carlo algorithm. 



From Eq. ( p.2| ) we can point out the natural dimensionless quantities to be used in the plot of 
our results: 

y = — n - and x = —=. (4.1) 

In Fig. [T] we plot our data all together, regardless of possible finite volume effects that are nev- 
ertheless smaller than statistical errors. We see that the continuum scaling law is satisfied for 
< x < 0.04. This allows us to extract the chiral condensate in the chiral limit, x — > 0, by the 
following linear relation: 

y = ao + a\x . (4.2) 



The best fit of the data to Eq. (4.2), in the scaling region, gives 



a = 0.07668 ± 0.00930 , a x = 11.20 ± 0.48 . (4.3) 

As a consequence, we see that in the chiral limit the external magnetic field does induce a non-zero 
chiral condensate. From Eq. ( |4.1[ ) and the result in Eq. (O), restoring cgs units, we find: 



— hceH 

(W) = — — (0.07668 ± 0.00930) . (4.4) 

271 

The non-zero value of the chiral condensate can be interpreted as the generation of a dynamical 



fermion mass md yn with gap Ao = md yn c . To extract Ao we use Eq. ( p.2| ) combined with Eq. (|4.4j): 



A ° - " V ^ (0-07668 ± 0.00930) . (4.5) 

iV/C(^)V 2tc 

5. Application to graphene 

The graphene is a single sheet of carbon atoms characterised by a honeycomb lattice described 
by two interpenetrating triangular Bravais lattices. The unit cell is instead a rhombus containing 
two carbon atoms. The reciprocal lattice comes out to be a hexagonal as well but, between the six 
corners of the Brillouin zone, only two points are inequivalent {i.e. not connected by a reciprocal 
lattice vector) and are labelled by K and K', called "valleys". 

A simple tight-binding model [ ]T2| ] shows that the dispersion relation E(k) goes to zero at K and 
K'. Everywhere else in momentum space E(k) / and we have two graphene bands, conduction 



4 



Lattice Planar QED in external magnetic field 



Pietro Giudice 



2.0 



1.5 - 



o 


L 




12 


ext 




1 


□ 


L 




12 


ext 




2 


o 


L 




16 


n t 

ext 




2 


A 


L 




16 


ext 




3 


V 


L 




24 


ext 




2 


< 


L 




24 


ext 




3 




± 



® 



V 



0.14 0.16 0.18 



0.20 



m Q / (eH/27c) 



Figure 1: Scaled chiral condensate versus the scaling variable x. The continuum line is the linear fit of the 
data in the scaling region < x < 0.04. 



and valence, lying symmetrically above and below the Fermi energy at E = 0. Therefore, the Fermi 
level is placed between the two symmetrical bands, with zero excitation energy needed to excite an 
electron from just below the Fermi energy to just above the six points of the corners of the Brillouin 
zone. 

Moreover, the dispersion relation turns out to be linear near the K and K' points: E(k) = 



±Hvf\k\, where vp is the Fermi velocity. As discussed in Ref. [|13|], this yields an effective Hamil- 
tonian made of two copies of Pauli spinors which satisfy the massless two-dimensional Dirac equa- 
tion with the speed of light replaced by the Fermi velocity. 

When graphene is immersed in a transverse magnetic field H the energy levels are quantised 
into non-equidistant Landau levels: 



v 



E n = sign(n)y2heH\n\-^, n = 0,±l,... . (5.1) 

The presence of the "anomalous" Landau level at zero energy leads to half-integer Quantum Hall 
(QH) effect; the Hall resistance is given by R xy = h/ (e 2 v), where the quantised filling factor v is: 

V = ±g,(AT+~) = ±2,±6,±10,... . (5.2) 
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The factor g s = 4 takes into account the spin and valley degeneracy and ± stands for electron and 
holes, respectively. 



Recently, it has been shown [ |l4] ] that for very strong magnetic field (H > 20 T) new QH states 
appear corresponding to v = 0, ±1 , ±4. This means that the n = Landau level is totally resolved 
into v = 0, ±1 plateaus, while the fourfold degeneracy in the n = ±1 Landau levels is only partially 
resolved into V = ±4 leaving a twofold degeneracy [jjjp. 

The new plateaus at v = 0,±4 can be explained by Zeeman spin splitting; v = ±1 can be 
explained if there is some kind of generation of a gap Ao, i.e. a valley symmetry breaking in n = 



Landau level [14, 15] 



In Ref. Jl5| ] the authors determine experimentally the AE(y = 1) dependence from the mag- 
netic field: they find for it a \[H behaviour. 



We have fitted the experimental data from Ref. Q15Q to 



AE(v = 1) = 2 (A (ff ) - |/i B #) , (5.3) 

where jXs is the Bohr magneton, g = 2 and Aq(H) ~ yH (note that the second term takes into 
account the Zeeman effect); we find: 



A (H) = (13.57 ±0.28) Kxk B yjH(T) , (5.4) 

where kg is the Boltzmann factor and H(T) stands for the magnetic field expressed in Tesla. 

It is believed that the generation of the gap is driven by the electron-electron interaction (for a 



review see Ref. [16]); in this picture, the value of the gap is expected to be of the order: 




Ao(H) « ~ 163Kxk B ^H(T) , (5.5) 



i.e. this result is one order of magnitude bigger than the experimental result of Eq. (p.4|). 

On the other hand, in Ref. [ [Tl| ] it is supported that the origin of the gap is related to a dynamical 
generation related to the presence of the magnetic field; this means that we can derive the value of 
the gap by our numerical results. 

The implied hypothesis here is that the Coulomb interaction between the electrons can be 
neglected because the dominant effect is the magnetic catalysis; this is the reason why we can apply 
planar QED directly to describe the graphene, although in the literature the Coulomb interaction is 
always considered as a 3d field acting on 2d fermions. 

To apply our results to graphene, we must simply replace the speed of light c with the Fermi 
velocity v F in Eq. (fO) [O]: 




A ~ - — ^— (0.07668 ± 0.00930) . (5.6) 

It is worth mentioning that the small gap approximation, which we are using here, is correct 
when it is satisfied the following relation: 



/ A ° < J2hv 2 F e/c rj A20K x k B ; (5.7) 
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the experimental value found in Eq. (5.4) is actually small compared with that in Eq. (5.7), therefore 
we can consider the small gap approximation valid in this context. 

Finally, using the experimental value for the Fermi velocity (vp « 1.0 x 10 8 cm/s) in Eq. (5.6), 
we get : 

A Q (H) = (5.52 ± 0.67) Kxk B \/H(T) . (5.8) 



This result can be compared directly with Eq. <J5T4|)- Since the two values are of the same order of 
magnitude, our result also confirm our hypothesis that the Coulomb interaction between electrons 
is actually negligible in presence of an external magnetic field. 

6. Conclusion 

Our nonperturbative Monte Carlo simulations have shown that the external magnetic field 
gives rise to a spontaneous breaking of the chiral symmetry. We think that it is possible to apply 
our results to describe the breaking of the valley symmetry in graphene under strong magnetic field 
and we determine the numerical value of the gap that compares quite well with the experimental 
value. 
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